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Abstract 

An interphase boundary may be immobilized due to nonlinear diffractional 
interactions in a feedback optical device. This effect reminds of the Turing 
mechanism, with the optical field playing the role of a diffusive inhibitor. Two 
examples of pattern formation are considered in detail: arrays of kinks in Id, 
and solitary spots in 2d. In both cases, a large number of equilibrium solutions 
is possible due to the oscillatory character of diffractional interaction. 



1 



Typeset using REVTJ5X 



There is a growing interest in transverse effects in nonlinear optics that manifest them- 
selves in spontaneous pattern formation in both active optical devices (lasers) and feedback 
systems (ring and Fabry - Perot cavities) driven by coherent sources. Symmetry breaking 
phenomena in passive nonlinear optical systems were studied in connection to self-focusing 
effects, stimulated Raman and Brillouin scattering and optical bistability in feedback sys- 
tems A balance of non-linearity and diffraction that is needed for pattern formation 
arises through sequential operation of the two effects in an optically thin nonlinear medium 
and in the empty part of the cavity. 

Transverse optical patterns have attracted particular attention in view of a possibility 
to imitate a variety of other pattern formation processes in non-equilibrium systems that, 
typically, require far longer observation times. Optical experiments are, however, at dis- 
advantage in that they have relatively small aspect ratios. The limiting factor here is the 
diffractional length that fixes the smallest scale of inhomogeneities of the intensity of the 
optical field. At the same time, far steeper inhomogeneities, limited by a short diffusional 
scale, are possible in a nonlinear optical medium. While the optical field cannot cause 
such inhomogeneities, it may stabilize them if they arise spontaneously, for example, near a 
Maxwell point in either equilibrium or non-equilibrium phase transition. 

The situation we envisage reminds of a common mechanism of formation of Turing 
patterns a nonlinear medium that can exist in two alternative states plays the role 

of a "slowly diffusing activator", while the role of the optical field is that of a "rapidly 
diffusing inhibitor". The analogy is, of course, not exact: nonlinear diffractive interaction 
in a feedback optical device is far more complex than diffusive interaction, and leads to a 
great number of alternative stable configurations. 

Our basic model is an optical ring cavity containing a thin sample of nonlinear medium 
with the thickness I in the direction of beam propagation negligible compared to the length L 
of the optical path in the empty part of the cavity. Under these conditions, diffraction of the 
light beam inside the nonlinear medium can be neglected, and consecutive transformations 
of the signal during a round-trip in the cavity include diffraction in the empty part of the 
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cavity, mixing with the pump beam, and point transformation in the nonlinear medium. 

In the paraxial approximation, the complex envelope amplitude E of the electric field in 
the empty part of the cavity obeys the parabolic equation 

iE z = V 2 E. (1) 

Here the coordinate z in the direction of propagation is scaled by the length L, and the 



transverse coordinates, by the diffraction length ^jL/k; V 2 denotes the transverse Laplacian. 

Since the material response time is typically much larger than the round-trip time in 
the cavity, we shall be interested in stationary solutions of Eq. ([I]) that give the field U(r) 
in the nonlinear medium corresponding to an instantaneous distribution of the material 
variable x(r). We define U as the field obtained by mixing the attenuated cavity signal 
Ef(r) = E(r, 1) with the field P of the pump beam. Then U(r) = P + aEf(r) where a is 
the transmission coefficient. The field E{ = E(r, 0) at the start of the next round trip is 
obtained by adding a phase shift /3% due to interaction with the medium. Expressing the 
solution of Eq. (JTj) by a functional Ef(r) = \& [-E'j(r)], and adding also a constant phase shift 
Q yields the functional equation for the quasistationary field U(r): 

U(r) = P + ae' n $ [u(r)e iMr) ] . (2) 

We adopt a model of material dynamics allowing for a phase transition biased by the 
intensity of the electric field in the medium. The simplest model of this type can be written 
as an evolution equation for the material variable x with a cubic nonlinearity, including a 
weak bias dependent on the intensity of the electric field: 

Xt = cr 2 V 2 X + x-X 3 + w(\U\)- (3) 

In the following, we shall use the bias function ty(|[/|) = {\U\ 2 — \Uq\ 2 ). The variable % can 
be interpreted as a phase field, in the spirit of common phenomenological phase-transition 
models || . The material response time is taken as unity; a is the ratio of the diffusional and 
diffractional lengths. Typically, a <C 1, so that the interphase boundary is sharp. As usual, 



we assume a thin sample approximation; then Eq.(|3]) retains only the transverse Laplacian 



V 2 . 



Eq. (0) is a dissipative dynamic equation of the form xt = —dF/Sx derivable from a 
suitable expression for the free energy T (playing the role of a Lyapunov functional): 



The last term, responsible for the field-induced bias, is recognized as a standard correction 
to the free energy of a dielectric in the electric field when the dielectric permittivity is 
linearly dependent on x- The model is therefore applicable wherever two dielectric phases 
with different permittivities may coexist. Experimentally, this situation can be realized, for 
example, near an isotropic - nematic transition. In this case, imposing a suitable orientation 
of the nematic phase can be used to enhance the difference between the refraction indices 
of the two phases. A more exotic possibility is a non- equilibrium phase transition between 
alternative steady states of a photochemical reaction. 

If U is eliminated with the help of Eq. (Q), slow transverse dynamics is determined by 
the evolution equation (^). The situation is somewhat similar to typical models of large 
aspect ratio dissipative pattern-forming systems where it is often possible to eliminate the 
"vertical" structure (in the direction of a flux sustaining the system in a non-equilibrium 
state). A substantial difference is that the effect of diffraction hidden in Eq. (0) is non-local, 
which makes the analysis far more difficult. 

Once the diffractional problem is solved, Eq. (|3]) can be treated in a standard way J7[^| . 
At a <C 1, X approaches one of alternative equilibrium values everywhere except boundary 
regions of 0(a) thickness (kinks). The phase equilibrium at \U\ = \Uq\ is biased by variable 
local field intensity, which sets the kinks into motion. The dynamics of a weakly bent kink 
is completely determined by its local geometry and the intensity of the biasing field \U\ at 
the kink position. The general equation of motion for a weakly curved kink is derived in a 
most transparent way with the help of multiscale expansion in a coordinate frame aligned 
with the kink ||. The kink propagates with the speed 
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c = —an =F bw(\U\), 



(4) 



where k is the curvature of the kink, and the coefficient b equals, for the cubic model in 
Eq.(JJ), to 3/y/2. The upper sign in Eq.(f|) applies to a kink with x(±oo) = ±1, and the 
lower, to an antikink with the reverse orientation. Curvature effects may balance the biasing 
field at k ~ w/a. This balance can be achieved when a curvature radius is measurable on a 
long ( diffract ional) scale, provided the bias is weak. 

In the following, I shall restrict to the dissipative limit when the transmission coefficient 
is small; the pump field is assumed to be constant in space. Then Eq. (|2]) yields, to the 
leading order in a, 



I shall consider first a one-dimensional model, and assume that the two phases coexist on 
an infinite x axis, with an interphase boundary (kink) at x — 0. In the zero order, the 
coexistence condition is \P\ = \Uq\. Assuming it holds, the profile of the material variable 
valid on the long (diffract ional) scale is simply xi x ) = 2H(x) — 1, where H(x) is the unit 
step function: H(x) = 1 at x > 0, H{x) = at x < 0. The jump is effected in a narrow 
0(a) interval that is negligible on the diffractional scale. Then Eq. (|1]) has to be solved with 
the initial condition 



The solution is expressed through the error function of a complex argument. This gives the 
cavity field u = u\(x) in Eq. (||) 



U(r) =P(l + cra(r)) 



(5) 



E(x,0)=P e" i/3 + 2i sin (3 H(x) 



(6) 




(7) 



Since Eq. (Q) is linear, solutions corresponding to n alternating kinks and antikinks are 



obtained by superposition: 



n 



i=i 



(8) 
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The corresponding bias of the field intensity is w = 2a\P\ 2 Keu n (x) . An 0(a) bias of the 
electric field at the kink location x = x k causes it to migrate, according to Eq. (f|), with the 
speed 

c = -3V2q k a\P\ 2 Reu n (x k ) (9) 

Incorporating various constants into the time scale brings the equation of motion of a 
system of kinks and antikinks to the form 



dx k 



qjqkip(x k - xj), ip(x) = Rewi(x). (10) 



This equation can be written in the gradient form 

dxk dV 



dt dx k jk 

where <§>(x) = \ ip(£)d£ 
Jo 



> V = ^q j q k <&(x k - Xj), (11) 



2 (x 2 + ix \ 
xi/j(x) + sin /3 sin Q) . (12) 

'7T V 4 / 



The function ip(x) is a combination of Fresnel integrals dependent on the phase angles 
(3 and VL. Both real and imaginary parts of the function erf(~y/ix) oscillate with increasing 
frequency and decreasing amplitude at \x\ — > oo around the asymptotic real value erf(oo) = 
1. These oscillations lead to a high complexity of solutions. Unlike long-range diffusive 
interaction in Turing patterns, the sign of diffractional interaction between kinks alternates 
with varying separation, and the potential Eq. ([12]) may have a great number of minima 
(see Fig. 1). 

We shall be mostly interested in the case when a kink and an antikink repel each other 
at short distances, which would lead to nucleation of kink-antikink pairs. This requires the 
derivative of the biasing field ^'(0) to be positive. Assuming < (3 < tt, this requires — j7r < 
Q < §7r. The repelling action may prevail, of course, on the diffractional scale only. Short- 
range attractive interaction takes over at 0(a) distances. The balance of attraction and 
repulsion determines the critical size of a nucleus. In one dimension, short-range interaction 
falls off exponentially ||10|| , and the size of the critical nucleus is a cr = 0(a In a) is very small. 



6 



Equilibrium separations of a kink/antikink pair, or, equivalently, the size of a solitary 
structure formed by the bound pair, are defined as zeroes of the biasing field ip(x). In a 
special case Q + [3 = ±7r/2, the asymptotic value of ip(x) at x — > oo vanishes. In this 
case, the equation ip(x) = has an infinite number of solutions. Asymptotically at large 
separations, 



By this formula, equilibrium separation distances are Xk = y (4/c — l)ir + Q, n — 1, 2, 

Alternating points correspond to stable equilibria. The asymptotic values are very close to 
"exact" solutions, except for xi, e.g. at Q — 0, f3 — ir/2 the exact value is x\ = 3.17, as 



Although the average interaction strength decays with separation not too fast (oc x^ 1 ), 
equilibrium positions of neighboring kinks in a multikink array are shifted only by few 
percentage points from positions computed taking into account pair interactions only. The 
minimal separation value for a symmetric triplet at Q = 0, /3 — vr/2 is ~ 3.35. There is, of 
course, also a variety of asymmetric equilibria forming an infinite grid at fl + f3 = ir/2. 

In multikink arrays the minimal separation between neighboring kinks tends to increase 
with the increasing size of the array. Early stages of evolution to equilibrium may appear 
to be rather disordered. An example of the relaxation process for a symmetric array of 11 
kinks, starting from positions corresponding to the minimal stable separation for a 3-kink 
array, is shown in Fig. 2. When the number of kinks is large, there are alternative close lying 
equilibria even when separations are close to minimal one. Such equilibria were detected 
when the dynamic equations of a symmetric 15-kink array were integrated starting from 
slightly different initial conditions. Multikink arrays may form chaotic structures containing 
an arbitrary number of gaps of different size. 

At fl + (3 differing from ±7r/2, the number of equilibria for kink-antikink pairs is finite, 
and is decreasing rapidly with growing deviation from the straight angle. The asymptotic 
estimate for a most distant equilibrium position is x max x (2/y / 7r)| sin/5/ cos(f2 + @)\. With 




(13) 
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no average bias \P\ 2 — \Uq\ 2 , equilibrium positions exist only while (3 remains within the 
range of maximal variation of this curve; e.g. at Q = there are no solutions at /3 > 0.657T 
or (3 < 0.417T. 

In two dimensions, an additional factor is the line tension, that affects the kink motion 
through the curvature term in Eq. (fj) . An island of one state (say, x ~ — 1) within a two- 
dimensional continuum of the alternative state would shrink, save for the optical interactions, 
when the phase equilibrium is unbiased. Due to nonlinear diffractional effects, the curved 
interphase boundary may stabilize at certain positions while retaining its circular shape 
because of the line tension. An a priori estimate for the radius a of a stationary island is 
a ~ a/a. 

Solving Eq. ([!]) in polar coordinates with the initial condition 



E(x, 0) = P e~ i/3 + 2i sin (3H{r - a) 



(14) 



yields the cavity field 



u r 



vyr, a 



~ ip + 2i sm(3v{r;a) 



1 — a I e 
h 



Ji(\a)Jo(\r)d\, 



(15) 



where J n (x) is a Bessel function. Equilibrium radii a verify the equation a 1 = QKeu(a), 
where Q = 2ba\P\ 2 a- 1 . At r = a, the integral in Eq. (|15| ) can be computed analytically: 

v(a; a) = l + J f y j exp f ^- j . (16) 



The resulting evolution equation of a is 

1 da Qh(a) sin/3 — 1 



a dt 
h{a) 



a 

,2 



a sin f h Q ) J ( — 

V 2 J \ 2 



Q cosVtcos (3\ 

2\ 



(17) 



A convenient form of the equation defining equilibrium values of a is 



{Qsm(3) 1 + a cos Vt cot (3 = h(a) 



(18) 



Solutions of this equation can be obtained as intersections of a straight line presenting its 
left-hand side with the plot of the function h(a) drawn in Fig. 3. The latter oscillates 
with an asymptotically constant amplitude between the maximum and minimum values 
^ cos 2 (|) ~ 0.963 and — sin 2 (|) « —0.165. The number of solutions is infinite when 
either (3 or Q equals to ±~, and the inverse of Q sin f3 lies within this interval. The stability 
condition with respect to radially symmetric perturbations is ah' (a) < h(a) — 1. At a 3> 1, 
stable solutions lie on descending segments in Fig. 3, i.e. at (2n — |)7r < a 2 < {2n — ■£). 

Stability against asymmetric perturbations dependent on the polar angle is tested || 
by perturbing the interphase boundary, r = a[l +C(<f>, t)], while retaining the quasistationary 
approximation for the optical field. If ((4>,t) ^ 1; the curvature of the kink is expressed as 
K = a~ 1 (l — ( — Cw)- Expanding the long-range biasing field w(r, cf>) in the vicinity of the 
unperturbed position brings Eq.({|) to the form 

dC/dt = a~ 2 a (Ch> + + b M", 4>) + w' s (a>)(] ■ (19) 

Here w' s (a) is the derivative of the stationary biasing field at the kink position obtained from 
Eq. flT5|), and w(r,<p) is the perturbation of the biasing field caused by the perturbation of 
the boundary. The latter field is obtained by solving Eq. with the initial condition 
E(r, (ft, 0) = 2i sin (3 ((4>)a(a). Setting £(0, t) = p n {t) cosn</> and adding, as usual, a constant 
phase shift, we obtain the scaled perturbation of the cavity field 

u(r, a) = liapn sin (3 e* n cosn0 « n (r, ; a); 

u n (r-a)= r Xe~ ix2 J n (\a)J n (\r)dX. (20) 
J o 

At n = 1, this expression is just the derivative du/dx of the stationary field Eq. fllS]) 
taken with the opposite sign; consequently, the two bracketed terms in Eq. fll9|) cancel, 
as they must do since this perturbation mode is equivalent to merely shifting the spot 
along the x axis. Computing the integral in Eq. ( p^| ) to express the perturbed biasing field 
w(a, 0) = 2a|P| 2 ReM(r, 0; a), and using it in Eq. ( |T9D yields the stability condition 

a~ 2 (n 2 - 1) > 2aQsm(3 [h n (a) - h^a)}, (21) 
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where 

POO 

hJa) = / A sin(A 2 - fi) J 2 (\a)d\ 
Jo 




At intermediate values of a, the quadrupole mode n = 2, deforming the disk into an ellipse, 
causes instability on the "stable" descending branches in Fig. 3 below the inflection point. 
At very large values of a, the solution is unstable against either n = 3 or n = 5 mode, 
depending on the value of Q. Instability may lead to the phenomenon of "replicating spots" 
similar to that observed in simulations of Turing patterns |Tl|. Spots may be arranged in 
crystalline arrays, though their shape remains nearly circular only at large separations. 
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FIG. 1. The relief of the potential for a kink - antikink - king triplet. 




FIG. 2. Evolution of separation distances between kinks and antikinks in a symmetric 11-kink 
array ((3 = 0, = 0). Darker curves correspond to kinks lying closer to the center. 
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